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Abstract 

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an arbiter concept 
considered for a possible Mars sample return mission. The primary experimental objectives were to 
characterize hypersonic near wake closure and determine if shear layer impingement would occur on 
the proposed arbiter afterbody at incidence angles necessary for a Martian aerocapture maneuvei . 
Global heat transfer mappings , surface streamline patterns , and shock shapes were obtained in the 
NASA Langley 20-Inch Mach 6 Air ' and CF 4 Tunnels for post-normal shock Reynolds numbers ( based 
on forebody diameter) ranging from 1 A00 to 415,000 , angles of attack ranging from -5 to 10 
degrees at 0, 3 , and 6 degree sideslip , and normal-shock density ratios of 5 and 12. Laminar , 
transitional , and turbulent shear layer impingement on the cylindrical afterbody was inferred from 
the measurements and resulted in a localized heating maximum that ranged from 40 to 75 percent of 
the reference forebody stagnation point heating. Comparison of laminar heating prediction to 
experimental measurement along the arbiter afterbody highlight grid alignment challenges 
associated with numerical simulation of three-dimensional separated wake flows. 


Nomenclature 

D aerobrake base diameter (in) 

d diameter (in) 

h heat transfer coeff. (lbm/ft 2 -sec), q/{ H aw - H w ) 

where H. IW = H t 2 
H enthalpy (BTU/lbm) 

M Mach number 

P pressure (psia) 

q heat transfer rate (BTU/ft 2 -sec) 

R radius (in) 

Re unit Reynolds number ( 1/ft) 

t time (sec) 

T temperature (°F) 

x axial distance along cylinder afterbody (in) 
a angle of attack (deg) 

<\> cone half angle (deg) 

0 shear layer turning angle (deg) 

p density (lbm/in 3 ) 

y ratio of specific heats 


Aerospace Technologist. NASA Langley Research Center, 
Hampton, VA. 

w Cooperative education student, NASA Langley Research Center. 
Hampton, VA. 

Copyright €'2001 by the American Institute of Aeronautics and 
Astronautics, Inc. No copyright is asserted in the United States 
under Title 17. U.S. Code. - The U.S. Government has a royalty-free 
license to exercise all rights under the copyright claimed herein tor 
government purposes. All other rights are reserved by the 
copyright owner. 


Subscripts 

aft vehicle afterbody 

eff angular dimension relative to a vector normal to 
base plane of aerobrake 
f forebody 

1 local condition 

oo free-stream conditions 

n model nose 

ref theoretical reference stagnation point heating 

s surface quantity, support sting 

t, 1 reservoir conditions 

2 stagnation conditions behind normal shock 

w wall 

Introduction 

As the result of a restructured 1 3 strategy for 
the exploration of Mars, it is anticipated that future 
robotic missions to Mars will continue to gather 
scientific data and demonstrate technologies required to 
address the feasibility of eventually establishing a 
human presence on Mars. The revised exploration 
architecture calls for a sample return opportunity as 
early as 201 1 4 . A Mars sample return (MSR) mission 
would collect terrestrial material from the planet surface 
and return these samples to Earth for detailed analysis. 
Although missions designed to return samples from 
Mars have been considered in the past 5 ' 7 , all planetary 
environmental and surface sample data to date has been 
acquired from remote based measurements obtained 
from the surface or in orbit of Mars. An MSR round- 


trip mission could require technologies such as 
aerocapture, precision landing, remote-controlled surface 
operations including surface material collection, orbital 
rendezvous, docking and sample transfer (all at Mars) 
culminating with a direct entry return at Earth (or 
shuttle rendezvous). 

Prior to the revised Mars exploration 
architecture, one approach considered for sample return 
called for a joint US/European MSR mission* 11 . Key 
elements of this mission profile are highlighted in Fig. 
1 . Collected surface samples would be placed into 
Mars orbit in advance of the arrival of a second vehicle, 
a French designed MSR orb iter (MSRO) : . Upon 
planetary approach, the MSRO vehicle would perform 
an aerocapture maneuver in the Martian atmosphere. 
With aerocapture. deceleration and orbital insertion is 
accomplished aerodynamically (in contrast to 
traditional propulsive techniques) using a blunt 
aeroshell. After an orbital rendezvous, the material 
samples would be autonomously transferred to the 
MSRO and stored in the US-designed Earth Entry 
Vehicle* (EEV) for the return transit to earth. 

The EEV and the retrieval hardware on the 
MSRO would reside behind the aeroshell during the 
aerocapture maneuver prior to its jettison. Proper 
positioning of this hardware is essential to insure 
aerodynamic stability of the MSRO during aerocapture, 
to reduce the chances of biological contamination of the 
EEV(s), and to avoid thermal damage from localized 
near wake phenomenon. Although it is generally 
recognized that convective heating rates to payloads 
behind aerobrakes are low 13 , localized heating 
maximum can occur if the boundary layer that separates 
from the shoulder of the aeroshell impinges on the 
afterbody 14 ' 17 . The complexity of hypersonic blunt 
body base flow is shown schematically in Fig, 2. 

The purpose of the present study was to 
experimentally investigate the near wake closure 
characteristics of a proposed MSRO vehicle. The 
primary objective was to experimentally determine if 
flow impingement occurred on the MSRO afterbody, 
and if present, determine the location and magnitude of 
the heating peak for comparison of prediction to 
measurement. Prediction methods that accurately 
model the near wake characteristics of blunt bodies are 
desired to reduce conservative design margins for 
planetary entry vehicles. 

Initial tests in the Langley 20-Inch Mach 6 
CF 4 Tunnel w ; ere used to assess laminar MSRO wake 
closure characteristics because of the facility's ability to 
provide normal shock density ratios (p 2 /p c< ,) of the same 
magnitude as that incurred in flight. It is anticipated 
that entry of the MSRO into the atmosphere of Mars 
(continuum-flow regime) during the aerocapture 
maneuver will produce values of shock density ratio 
near 20 (these high values of normal shock density- 
ratio encountered in hypervelocity flight are produced 


from dissociation of atmospheric gases as they cross 
through the shock wave into the shock layer and are 
often referred to as real gas effects). Flight values of 
normal shock density ratio are significantly larger (3 to 
4 times) than those produced in conventional blowdown 
hypersonic wind tunnels using air or nitrogen as a test 
gas lx . The NASA Langley 20-Inch Mach 6 CF 4 
Tunnel utilizes a test gas (tetrafluoromethane-CF 4 ) with 
a molecular weight 3 times that of air in order to 
generate a normal shock density ratio of 12 thereby 
simulating this particular aspect of a real gas. The 
simulation is achieved at moderate levels of enthalpy 
and without dissociation of the test gas. It is well 
recognized the normal shock density ratio is one of the 
primary flight simulation parameters that govern the 
inviscid flow and aerodynamics of blunt bodies at 
hypersonic speeds 1 ^. Complimentary tests in a perfect 
gas Mach 6 air facility were later performed over a 
larger Reynolds number range to assess the effects of 
shear layer transition on wake closure and impingement 
heating. 

The MSRO was tested at normal shock 
density ratios of 12 and 5 (in CF 4 and air. respectively) 
over a range of freestream unit Reynolds number from 
1 2 x 10 4 to 8 x 10 6 per foot (post shock Reynolds 
number range based on aeroshell diameter of 1,400 to 
415,000). Angle of attack was varied from -5 deg to 10 
deg at zero deg sideslip (limited data at 3 and 6 deg 
sideslip). Test techniques that w'ere utilized during 
these tests include thermographic phosphors, which 
provided global quantitative surface heating images; oil- 
flow, which provided surface streamline information; 
and schlieren, which provided shock system details. 
Emphasis was placed on the afterbody surface heating 
augmentation due to shear layer impingement. 

Experimental Methods 

Model 

The Outer Mold Lines (OML) of the proposed 
MSRO vehicle shown in Fig. 3 represents a reference 
baseline concept for which experimental and 
computational wake assessment studies can be 
conducted. The proposed MSRO vehicle consists of a 
drag brake and a base-mounted cylindrical afterbody 
(d att /d t =0.48). The base of the drag brake (aeroshell) is 
concave inwards as shown in Fig. 3. A single EEV for 
sample storage and a housing structure for the 
retrieval/transfer hardware are positioned on the 
afterbody. The drag brake, or aeroshell, is based upon a 
design originally proposed in support of the Aeroassist 
Flight Experiment (AFE) :o - a 60 degree elliptically 
blunted cone with the base raked at a 73 degree angle. 
At zero degrees angle of attack, the wind vector is 
aligned with the minor axis of the ellipsoidal nose (see 
Fig. 3). 

A photograph of tw'o 4-in. diameter (0.028 
scale) cast ceramic MSRO models used in the test is 
shown in Fig. 4. For construction of the fused silica 
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models, a rapid prototyping technique was utilized. 
Stereolithography (SLA) was used to build a resin 
aeroshell model with various, detachable SLA 
afterbodies representing a component buildup of the 
return orbiter EEV and retrieval/transfer housing 
structure. The SLA resin aeroshell was then assembled 
with the desired afterbody; the mated pair served as a 
pattern to construct molds from which the cast ceramic 
MSRO models were made. A support sting enters 
from the rear of the cylindrical afterbody. A magnesia 
ceramic was used to backfill the cylindrical afterbody, 
thus providing strength and support to the sting 
support structure. The aeroshell baseplanes on the 
MSRO models were constructed both open (concave 
inwards as sketched in Fig. 3) and closed. A 
photograph of these two base configurations is shown 
in Fig. 4. 

Typically, two casts of each configuration are 
made; the primary being immediately prepared for 
testing and the back-up shell held in reserve, in the 
event of problems with the primary. In order to obtain 
accurate heat transfer data with the phosphor technique, 
the models are cast with a material having low thermal 
diffusivity and well defined, uniform, isotropic thermal 
properties. The phosphor coatings typically do not 
require refurbishment between runs in the wind tunnel 
and have been measured to be approximately 0.001 
inches thick. Details concerning the ceramic model 
fabrication technique and phosphor coating can be found 
in Refs. 21 and 22. Fiducial marks were placed on the 
model surface to assist in determining spatial locations 
accurately. 

Once the phosphor testing was completed, the 
untested backup ceramic models were used for 
visualization tests (oil-flow and schlieren). 

Facility Descriptions 

Tests were conducted in two facilities managed 
under the Aerothermodynamic Facilities Complex 
(AFC). This complex presently consists of four 
hypersonic wind tunnels that represent a large fraction 
of the nation's conventional aerothermodynamic test 
capability 2 '. Collectively, they provide a wide range of 
Mach number, unit Reynolds number, and normal 
shock density ratio. This range of hypersonic 
simulation parameters is due, in part, to the use of two 
different test gases (air and tetraflouromethane), thereby 
making the facilities unique national assets. The AFC 
facilities are relatively small and economical to operate, 
hence ideally suited for fast-paced aerodynamic 
performance and aeroheating, and transition studies 
aimed at screening, assessing, optimizing, and bench- 
marking (when combined with computational fluid 
dynamics) advanced aerospace vehicle concepts and 
basic fundamental flow physics research. 

20-Inch Mach 6 Air Tunnel: Heated, dried, 
and filtered air is used as the test gas. Typical 
calibrated operating conditions for the tunnel are: 


stagnation pressures ranging from 30 to 500 psia; 
stagnation temperatures from 760 to 1000-degree R; 
and freestream unit Reynolds numbers from 0.5 to 8 
million per foot. A two-dimensional, contoured nozzle 
is used to provide nominal freestream Mach numbers 
from 5.8 to 6.1. The test section is 20.5 by 20 inches; 
the nozzle throat is 0.399 by 20.5-inch. A bottom- 
mounted model injection system can insert models 
from a sheltered position to the tunnel centerline in less 
than 0.5-sec. For the transient heat transfer tests, the 
model residence time in the flow' was limited to 5 
seconds. A detailed description of this facility may be 
found in Ref. 23. 

Mach 6 CFj Tunnel : Heated, dried, and 
filtered tetrafluoromethane (CF 4 ) is used as the test 
gas. Typical calibrated operating conditions for the 
tunnel are: stagnation pressures ranging from 85 to 
2000 psia, stagnation temperatures up to 1300-degree 
R,and freestream unit Reynolds numbers from 0.01 to 
0.41 million per foot. A contoured axisymmetric 
nozzle is used to provide nominal freestream Mach 
numbers from 5.9 to 6.2. The nozzle exit diameter is 
20 inches with the flow exhausting into an open jet 
test section; the nozzle throat diameter is 0.466-inch. 
A bottom-mounted model injection system can inject 
models from a sheltered position to the tunnel 
centerline in less than 0.5-sec. For the transient heat 
transfer tests, the model residence time in the flow' was 
limited to 5 seconds. A detailed description of this 
facility may be found in Ref. 19. 

Test Conditions and Setup 

Nominal reservoir and corresponding free 
stream flow conditions for the MSRO test series are 
presented in Table 1. The freestream properties w'ere 
determined from the measured reservoir pressure and 
temperature and the measured pitot pressure at the test 
section. Static pressure in the CF 4 test section that 
enclosed the open jet were monitored to assess the 
potential for contraction of the open-jet test core flow 
with time during any given run. The ratio of pitot to 
reservoir pressure was also measured in both facilities 
during each run to determine flow' conditions and test 
times. The ratio of projected model frontal area to core 
flow cross sectional area for the 0.028-scale model is 
approximately 0.15. 

All models were supported by a cylindrical steel 
rod (sting) which extended downstream from the model 
base at a 105-degree incidence angle relative to the 
forebody base plane. This was done to align the sting 
along a computationally predicted wake axis with the 
model at zero degree incidence. To determine the 
sensitivity of the support system on base flow heating, 
several runs were made with an adapter sleeve placed 
over the existing sting which increased the sting to 
forebody diameter ratio (d/dd from 0.19 to 0.32. 
Model angle-of-attack and sideslip were set to zero in 
the tunnel using a combination of an inclinometer and 
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a laser alignment system. A photograph of the 
installation in the CF 4 tunnel is shown in Fig. 5. 

Test Techniques 

Global Phosphor Thermography: Advances 
in image processing technology which have occurred in 
recent years have made digital optical measurement 
techniques practical in the wind tunnel. One such 
optical acquisition method is two-color relative- 
intensity phosphor thermography 24 ' 26 which has been 
utilized in several aeroheating tests conducted in the 
hypersonic wind tunnels of NASA Langley Research 
Center. 27 70 . With this technique, ceramic wind tunnel 
models are fabricated and coated with phosphors that 
fluoresce in two regions of the visible spectrum when 
illuminated with ultraviolet light. The fluorescence 
intensity is dependent upon the amount of incident 
ultraviolet light and the local surface temperature of the 
phosphors. By acquiring fluorescence intensity images 
with a color video camera of an illuminated phosphor 
model exposed to flow in a wind tunnel, surface 
temperature mappings can be calculated on the portions 
of the model that are in the field of view r of the camera. 
A temperature calibration of the system conducted prior 
to the study provides the look-up tables that are used to 
convert the ratio of the green and red intensity images 
to global temperature mappings. With temperature 
images acquired at different times in a wind tunnel run, 
global heat transfer images are computed assuming one- 
dimensional semi-infinite heat conduction. The 
primary advantage of the phosphor technique is the 
global resolution of the quantitative heat transfer data. 
Such data can be used to identify the healing footprint 
of complex, three-dimensional flow phenomena (e.g., 
separation/reattachment, transition fronts, turbulent 
wedges, boundary' layer vortices, etc.) that are extremely 
difficult to resolve by discrete measurement techniques. 
Because models are fabricated and instrumented more 
rapidly and economically, global phosphor 
thermography has largely replaced discrete heating 
instrumentation in Langley’s AFC. 

Schlieren photography: Flow visualization 
in the form of schlieren was used to complement the 
surface temperature and heating measurements. The 
LaRC 20-lnch Mach 6 Air and CF 4 Tunnels are 
equipped with a pulsed white-light, Z-pattern. single- 
pass schlieren system with a field of view 
encompassing the entire test core. The light sources 
are pulsed for approximately a 3 ms. Images were 
recorded on a high-resolution digital camera and 
enhanced with commercial software. 

Oil Flow Visualization: Surface streamline 
patterns were obtained using the oil-flow technique. 
Backup ceramic models were spray-painted black to 
enhance contrast with the white pigmented oils used to 
trace streamline movement. A thin basecoat of clear 
silicon oil was first applied to the surface, and then a 
mist of pinhead-sized pigmented-oil drops was applied 


onto the surface. After the model surface was prepared, 
the model was injected into the hypersonic flow, and 
the development of the surface streamlines was recorded 
with a conventional video camera. The model was 
retracted immediately following flow establishment 
(and formation of streamline patterns) and post-run 
digital photographs were taken. 

Data Reduction and Uncertainty 

A 16-bit analog-to-digital facility acquisition 
system was used to acquire flow' condition data. 
Measured values of P t j an( j T t j are believed to be 
accurate to within ±2 percent. 

Heating rates were calculated from the global 
surface temperature measurements using one- 
dimensional semi-infinite solid heat-conduction 
equations, (see Ref. 26). As discussed in Ref. 26, the 
accuracy of the phosphor system measurement is 
dependent on the temperature rise on the surface of the 
model. For the heating measurements, the phosphor 
system measurement accuracy is believed to be better 
than ±8%, and the overall experimental uncertainty of 
the phosphor heating data due to all factors is estimated 
to be ±15%. In areas on the model where the surface 
temperature rise is only a few degrees (i.e. cylinder 
afterbody, upstream of flow reattachment), the 
estimated overall uncertainty is on the order of ±25% . 
Repeatability for the normalized centerline (laminar) 
heat transfer measurements on the afterbody was found 
to be generally better than ±8%. Uncertainties in 
model angle-of-attack and sideslip are believed to be 
±0.2 degree. 

Computational Methods 

The Langley Aerothermodynamic Upwind 
Relaxation Algorithm (LAURA) 71 72 was used to 
provide laminar numerical simulations of the MSRO 
wake structure at CF 4 wind tunnel conditions. 

Grid topology, sensitivity studies, and convergence 
metrics are discussed in Ref. 1 1 . On the surface, no-slip 
conditions were applied and model wall temperature was 
assumed isothermal and set to 80 deg F (T w /T=0.43). 

Results and Discussion 

Preface 

The heating mappings and distributions in the 
subsequent sections are presented in the form of a 
normalized heating coefficient, h/h rc f where h rc { 
corresponds to the theoretical stagnation point heating 
at oc=0 degrees on the model forebody using the method 
of Ref .33. The reference radius is 1.81 -inches, which 
approximates the ellipsoidal nose radius of the AFE 
forebody in the model plane of symmetry. For the 
global images, the colors tending towards red indicate 
areas of higher heating (temperatures) while the colors 
towards blue represent regions of lower heating. Shear 
layer turning angle, which is measured relative to a 
direction normal to the base plane of the aerobrake (see 
Fig. 3), is inferred from schlieren images, surface 
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streamline patterns, and local heating. Unless indicated 
otherwise, all data presented was obtained on models 
with the aeroshell baseplane configured as open (see 
Fig. 4). 

Shear layer Identification 

The MSRO is intended to trim between a=-4 
and +4 degrees during its aerocapture maneuver at Mars. 
Schlieren images associated with the MSRO flow field 
(in M w =6 air at Re 2n =4.15 x 10 5 ) are presented in Fig. 
6, for a range of angle of attack (-5<ot<+10 degrees). 
In the sequence of images, a weak lip shock associated 
with the rapid over expansion and subsequent 
recompression around the aerobrake shoulder is evident, 
in addition to a relatively thin free shear layer. While 
not attempted in the present study, M^b air flow' field 
surveys in the base region of a similar blunt body 14 
have correlated wake lip shock and shear layer 
boundaries inferred from pitot pressure profiles 
to visual locations inferred from schlieren 
images. The MSRO schlieren images shown in 
Fig. 6, indicate that at angles of attack larger than +6 
degrees, the afterbody upper symmetry plane (surface 
opposite the EEV and retrieval hardware) avoids direct 
shear layer impingement for this condition. As the 
model angle of attack is decreased, the onset of shear 
layer impingement is observed to occur (in air) between 
ot=+6 and +5 degrees, and is indicated by a weak 
recompression shock at the reattachment point near the 
end of the cylindrical afterbody. The shear layer 
impingement point progressively moves towards the 
aerobrake base (and the strength of the recompression 
shock increases) as the vehicle angle of attack is 
decreased to a=-5 degrees. At the limiting incidence 
angle of a=-5 degrees, an interaction of the lip shock 
with the afterbody recompression shock is observed 
near the aft comer of the afterbody. Shear layer 
impingement on the EEV and retrieval hardware are not 
observed in M^=6 air for -5<a<+10 degrees. The 
sensitivity of the schlieren system is insufficient to 
detect the wake shear layer in the CF4 facility due to 
the low density of the base region flow . 

Heating distributions obtained in M w =6 air and 
CF 4 along the cylinder afterbody at several time 
intervals are used to determine when the base flow is 
established. Although some degree of wake flow- 
unsteadiness is inherent, the shear layer position is 
typically stable within the framing rate (30 frames/sec.) 
and spatial resolution of the schlieren video imagery. 

Shear layer Impingement 

Shear layer impingement will result in localized 
surface pressure and convective heating increases. 
Typical surface streamlines and heating distributions on 
the MSRO afterbody in the presence of shear layer 
impingement are shown in Figs. 7 and 8 at (X=-4 
degrees for M^=6 air and Re 20—6 .5 x 10 4 . Similar 
trends were observed at corresponding conditions in 


M^=6 CF 4 . The pressure rise due to reattachment and 
downstream recompression contributed to the reversal 
of the flow back towards the open forebody base as 
inferred from the surface streamline patterns (see Fig. 7). 
The location of the dividing streamline delineating the 
entrained recirculating base flow from the flow 
continuing downstream is evident in this figure. As 
discussed in a subsequent section, this reversed flow 
was often observed to produce complex recirculation 
patterns on the forebody base presumably due to 
primary and secondary vortices. The corresponding 
surface heating (see Fig .8), indicates that the local 
impingement heating at this condition approached 65% 
of the reference forebody stagnation point value; 
heating levels upstream of the attachment location and 
on the base were under 5% of the reference stagnation 
value. Regions of flow reattachment inferred from 
heating patterns on the afterbody side is not evident for 
sideslip angles up to 6 degrees (not shown). 

MSRO afterbody wall temperatures in M«=6 air 
obtained 5 seconds after model exposure were generally 
uniform (T w =l 10 deg. F; T u /T t =0.6n in locations away 
from impingement. At higher Reynolds numbers, 
localized turbulent heating near impingement produced 
wall temperatures as high as 230 degree. F (T w /T t =0.74) 
with peak heating levels that approached 75 % of the 
reference stagnation value. 

Reynolds Number Effects 

Because experimental evidence revealed shear 
layer impingement at trim angles of attack, it became 
necessary to quantify impingement heating levels and 
thus identify the state of the prevailing interaction 
(laminar, transitional, or turbulent). As discussed in 
detail in Ref. 34, separated flows are largely 
characterized by the prevailing boundary layer or shear 
layer state. Differences in laminar and turbulent 
separated regions are primarily attributed to enhanced 
mixing (greater effective momentum transfer) associated 
with turbulence. The presence of a transitional shear 
layer often leads to much higher peak heating at 
reattachment than if the separated flow is entirely 
laminar or turbulent 1435 . While edge conditions (M. 
Re, and y) at the point of boundary layer separation 
determine the initial shear layer turning angle, the state 
of the interaction (reattachment) downstream also 
effects the extent (and degree) of unsteadiness of the 
reversed flow. Thus, a secondary' objective ot the 
experimental program was to produce both a laminar 
and turbulent interaction from which comparisons to 
numerical solutions could be made. Confidence levels 
associated with numerical simulation of separated 
laminar flows are generally lower than for laminar 
attached flows. The additional complexity of accurately 
modeling transitional or turbulent separated flow has 
also proven to be a challenge 36 (ironically, it is the 
laminar wake interaction that is often difficult to 
achieve in conventional ground based hypersonic wind 
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tunnel due to facility operational limitations at low 
Reynolds numbers). 

The M & =6 CF 4 tunnel was initially attractive 
from the perspective of testing in a gas with a low ratio 
of specific heats (y) - as is encountered in 
hypervelocity flight. In addition, the facility provides 
the best opportunity from which to maintain a laminar 
interaction since it operates at low pressures (and hence, 
Reynolds numbers). Normalized heating along the 
MSRO cylindrical afterbody is shown in Fig. 9 at a=-4 
degrees for M„=6 CF 4 over a range of post normal 
shock Reynolds numbers (Re 2D =I ,420 to 52,442). The 
abscissa of Fig, 9 (and all subsequent afterbody heating 
distributions) refers to the running length along the 
cylinder non-dimensionalized by the forebody diameter, 
originating at the aft comer and running forward 
towards the aerobrake base. The relatively small 
invariance in the measured magnitude and spatial 
location of the reattachment heating peak for 
Re :n = 1,420 to 3.946 is interpreted as evidence of a 
laminar interaction. An increase (and subsequent 
movement) of peak reattachment heating towards the 
aerobrake base for Re 2D >3,946 suggested a non-laminar 
interaction 14163 ' and/or an increase in the shear layer 
turning angle. The effects of Reynolds number on 
shear layer turning angle will be discussed in a 
subsequent section. 

The M.=6 air tunnel provides a two order-of- 
magnitude increase in post- normal shock Reynolds 
number over the CF 4 data. Normalized heating along 
the MSRO cylindrical afterbody is shown in Fig. 10 at 
a=-4 degrees in M^=6 air for Re :D =36.700 to 415,000. 
A small overlap in terms of post-normal shock 
Reynolds number exists between the two facilities. 
Qualitatively, the magnitude of the measured 
impingement heating peak in M^=6 air at the lowest 
Reynolds number (Re 2D =36,700- Fig. 10) is consistent 
with that measured at the highest Reynolds number 
(Re 2D =52,442-Fig. 9) in M oc =6 CF 4 . The continuation 
of the forward movement of the local reattachment 
heating peak in M.^6 air toward the aerobrake base is 
evident. The continual rise in the heating peak to 
h/h ret -0.66, as measured at Re :r) =65.000 and the 
subsequent drop to h/h re p=0.60 for Re 2D =123,000, (see 
Fig. 10), is interpreted as the condition at which the 
transitional shear becomes fully turbulent. It is 
suggested that this situation is analogous to the 
transitional “overshoot” phenomenon commonly 
observed with attached wall boundary layers (where it 
has been conjectured that the larger vortical length 
scales associated wdth transitional flows are more 
effective at momentum transfer than the finer scales 
found in a fully turbulent situation). In the postulated 
turbulent regime, a linear rise in the non-dimensional 
heating peak to 0.73 for Re :D =4 15,000 w'as measured in 
M^=6 air. A Reynolds number collapse of turbulent 
heating data is not observed when normalized with a 
laminar reference value. 


The approximate range of Re :n at which the 
wake interaction has been postulated to go from a 
laminar to a transitional state, and from a transitional 
state to a fully turbulent are presented in Fig. 11. In 
this figure, the local heating peak at reattachment has 
been plotted as a function of Re :n for both M & =6 air 
and CF 4 . As inferred in Fig. 9. the relatively small 
invariance in the measured magnitude of the 
reattachment heating peak for Re 2D =l ,420 to 3,946 was 
inteipreted as evidence of a laminar interaction in M^=6 

cf 4 . 

A limited literature search yielded several 
empirically derived correlations from which to 
determine the onset of wake flow transition. In a 
qualitative attempt to classify the present M„=6 CF 4 
results as laminar, the data are presented in terms of. 
(and compared against), a simple unified wake 
transition correlation^. Fig. 12. The parameters for the 
original correlation in Ref. 38 are freestream length 
Reynolds number (the length based upon the axial wake 
transition location relative to the base) and freestream 
Mach number. The correlation was developed from a 
large experimental database for cones, spheres, wedges, 
and sphere-cones that were tested over a range of 
supersonic and hypersonic Mach numbers. The 
potential effects of bluntness and body shape on w^ake 
transition were recognized in Ref. 38 and the length 
Reynolds number parameter w^as modified to include the 
square of the ratio of the freestream Mach number to a 
local wake Mach number. In terms of the correlation, 
the present M^-6 CF 4 condition (interpreted as laminar 
for Re 2D =l ,420) is an order of magnitude or more below r 
the specified criteria for wake transition. The local 
wake Mach number for the present results was 
determined from laminar prediction corresponding to a 
location just outside of the recirculating base tlow (near 
shear layer impingement). The size of the symbols, 
(Fig. 12), are indicative of the dependence of the 
correlation parameter to the variation in Mach number 
as computed in this region. The transition length was 
assumed to be the axial distance of the aft cylinder 
(presumably where transition was observed) to the 
MSRO base. The M w =6 CF 4 results which were 
interpreted as transitional (Re 2D =52,442) were 
significantly closer to the correlation criteria for wake 
transition. 

Base Effects 

It is generally recognized that the presence of a 
support sting will have an effect on the base flow 
characteristics of a blunt body in hypersonic flow . A 
limited data set was taken to assess support interference 
effects. Experimentally, the sensitivity of the MSRO 
afterbody heating to the potential effects of a cylindrical 
model support sting was assessed for sting-to-model 
forebody diameter ratios (djd^) of 0,19 and 0.32. Fig. 
13 shows that at a transitional Reynolds number 
(Re 2D =36.700), for M^b air at a=-4 degrees, an 
increase in the support sting diameter produces no 
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measurable effect. Transitional data for lS/1^6 CF 4 (not 
shown) supports a similar conclusion. Comparison of 
the data in Fig. 13 also serves to indicate the level of 
repeatability associated with the wake heating 
measurements. 

As shown in Figs. 14, 15, and 16, the state of 
the near wake (laminar, transitional, and turbulent 
respectively) has a pronounced effect on the magnitude 
of the peak heating associated with shear layer 
impingement for MSRO models (both open and closed 
baseplanes-see Fig. 4). For a laminar wake condition 
(Re 2D =3,946 ) in M r .=6 CF 4 at a=-4 degrees, the closed 
(or filled baseplane) appeares to alter the base flow and 
lower the afterbody peak heating relative to that 
obtained with the open base by approximately 60% 
(from h/h re f=0.43 to 0.27), Fig. 14. In contrast, a 
turbulent wake condition (Re 2n =4 15.000) in 1^=6 air 
increases afterbody peak heating for the closed 
baseplane by approximately 23% (from h/h re p0.73 to 
0.90), Fig. 16. Evidence of the complexity of the 
turbulent reversed flow impingement onto the closed 
baseplane at M m =6 air (Re :n =4 15,000) is shown in the 
surface streamline patterns in Fig. 17. Baseplane 
surface heating in this vicinity (not shown) is 
approximately h/h rc j— 0. 1 2. When the shear layer state 

is interpreted as transitional in VL=6 air (Re 2D =52,442), 
no changes in peak heating magnitude or location for 
the open and closed baseplanes is observed. Fig. 15 (a 
similar conclusion was obtained for transitional 
measurements in M oc =6 CF 4 - not shown). Thus, the 
turbulent wake associated with the closed MSRO base 
represents a worst case scenario in terms of afterbody 
heating where localized values at reattachment 
approached forebody stagnation levels. 

Shear layer Turning Angle Characteristics 

From a practical perspective, the usable volume 
for payload placement behind an aerobrake is generally 
constrained by the position of the wake free shear layer. 
The wake boundaries from which the MSRO afterbody 
was initially designed w^ere inferred from a correlation 
derived from a series of ground-based blunt body tests at 
M.^6 and 10 air and numerical flight prediction 3 ’. In 
this correlation, measured and computed shear layer 
deflection angles and vehicle angle of attack are 
expressed relative to a direction normal to the given 
aerobrake baseplane. A linear relationship was 
identified between the shear layer turning angle (0 eff ) 
and the effective angle of attack (a ctY ). 

The shear layer turning angle (as a function of 
effective angle of attack) inferred from the present 
heating, surface streamline, and shock pattern 
measurements in M m =6 air and CF 4 are presented in 
Fig. 18 and Fig. 19, respectively. For clarification 
purposes, the ground based results from Ref. 39 are 
presented in these figures in the form of a shaded band. 
Consistent with the methodology of Ref. 39, the shear 
layer deflection angles are measured relative to a 


straight line drawn from the model base comer to the 
measured location of the attachment point as inferred 
from surface streamline patterns, schlieren images, and 
the location of peak heating using phosphor 
thermography. 

The present MSRO results obtained with the 
open base plane exhibit the same general linear trend as 
the correlation of Ref. 39. The turning angle interred 
from the CF 4 heating and streamline measurements (see 
Fig. 19), compare well with the correlation. In the 
present data, the variation of the turning angle for each 
effective angle of attack has been attributed to 
transitional/turbulent Reynolds number effects 
discussed earlier. That is, for a non-laminar interaction, 
an increase in Reynolds number moves the 
impingement location towards the aerobrake base, 
implying a larger turning angle. The effect of 
Reynolds number on shear layer turning tor a laminar 
interaction is inconclusive due to the limited range of 
Reynolds number for laminar CF 4 data. The M,.=6 air 
results. Fig. 18, indicate a bias towards a huger turning 
angle relative to the original correlation and this is. in 
part, attributed to the greater turbulent Reynolds 
number range achieved in the present tests. 

For a given Re^ n corresponding to a non-laminar 
shear layer, the shear layer turning angle as interred 
from surface streamlines or shock patterns is generally 
greater than that determined via the peak heating 
location. Shear layer turning angles as a function of 
Re :o for a=-4 degrees in M & ,=6 air and CF 4 are 
presented in Fig. 20 and Fig. 2 1 . respectively. For the 
laminar M w =6 CF 4 condition (Re 2D = 1,420 to 3,946), 
the turning angle associated with peak heating on the 
afterbody is within one degree of that inferred from 
surface streamlines or shock patterns, as shown in Fig. 
21. For transitional/turbulent conditions in M x =6 air 
(Re 2D > 3,946, Fig. 20), this difference is as large as 
four degrees. As indicated by the data in Figs. 20 and 
21, the turning angle for a transitional/turbulent 
interaction is larger than its laminar counterpart. As 
inferred from this data, the transitional/turbulent 
heating maximum on the cylinder afterbody 

(Re :D >3,946) is located downstream of the reattachment 
point as determined from surface streamlines or shock 
patterns. 

It was concluded in Ref. 39 that no significant 
effects due to M«, Re :[> or vehicle forebody geometry 
were observed in the ground-based test data that the 
shear layer turning angle correlation w r as derived from. 
Slightly larger shear layer turning angles were predicted 
at flight conditions and it was suggested that this could 
have been a gas chemistry effect. As discussed in Ref. 
39, the severe compression in the forebody shock layer 
and the subsequent rapid expansion into the w p ake are 
characterized by high temperature, thermo-chemical 
nonequilibrium processes. In contrast to quantifying 
the aerodynamic effects associated a real gas, the 
applicability of the present MSRO simulation tests in 
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CF 4 are not certain regarding the characterization 
ofblunt body wake flow. From an aerodynamic 
perspective, it has been shown 40 that the increments and 
trends provided by real gas simulation tests in air and 
CF 4 are applicable to flight provided that ( 1 ) the vehicle 
aerodynamics are dominated by the windward surface 
pressure (shear stress contributions negligible), (2) y 
within the flight windward shock layer does not 
significantly vary spatially, and (3) y within the flight 
windward shock layer is close in magnitude to that 
produced in CF 4 (y=l .1 ). Despite ground based testing 
simulation issues and geometrical afterbody differences 
from those configurations used to define the turning 
angle correlation of Ref. 39, the present MSRO 
laminar CF 4 data has suggested the general applicability 
of the correlation for predicting wake shear layer 
turning for initial design estimates. If shear layer 
transition in flight is anticipated (i.e. via forebody 
ablation or surface roughness) a larger turning angle 
should be expected. 

The intent of the baseline MSRO afterbody 
design was to avoid impingement at trim angles of 
attack. The present results along with recent 
hypersonic tests 1241 conducted in European facilities (in 
air and CO : ) have all indicated flow impingement on 
the MSRO afterbody near trim conditions. Flight 
computations 10 12 have indicated impingement as well. 
These experimental and computational results reveal the 
presence of impingement w hich w ill necessitate the use 
of a supplemental thermal protection system such as an 
afterbody shield discussed in Ref. 12. In hindsight, it 
appears that the early engineering estimates of the 
MSRO wake closure characteristics (to avoid shea:' 
layer impingement) were not conservative enough. For 
design purposes, a shear layer turning angle of 30 
degrees (see Fig. 19) was originally assumed based 
upon an effective angle of attack of 15 (oc=+2 degrees) 
degrees. As aerodynamic performance of the aerobrake 
during an aeropass maneuver was refined to include 
dispersions, the effective angle of attack to trim the 
MSRO was increased to 19 degrees (a=-2 degrees) and 
later to 21 degrees (a=-4 degrees). Based upon the 
correlation of Ref. 39. the projected shear layer turning 
angle relative to the baseplane would increase to 
approximately 39 degrees. The relationship of these 
turning angles relative to the MSRO afterbody is 
shown. Fig. 22. As inferred from the correlation, 
impingement on the afterbody would occur for effective 
angles of attack of 19 and 21 degrees. A small rotation 
or lateral displacement of the afterbody in the plane of 
symmetry (or scarfing the aft comer) would possibly 
provide a means to avoid impingement. An afterbody 
shield would also serve to prevent impingement but 
would add undesirable mass to the orbiter. 

Comparison to Measurement 

Accurate. quantitative measurement and 
numerical simulation of laminar, transitional, and 
turbulent separated/re attaching flows remains a 


challenge. Due to the operating range of most 
continuum hypersonic facilities, maintaining a laminar 
shear layer in regions of strong expansion (such as that 
encountered in the near wakes of blunt bodies) has in 
the past proven difficult 14 The influence of tunnel 
noise on shear layer transition is not well understood. 
When turbulent interactions are present, numerical 
turbulence models play a crucial role in the simulation 
of these non-laminar flow’s where separation, shock- 
boundary layer interaction, and flow reattachment are 
present 10 . For the simpler laminar interactions, w'hich 
avoid turbulence modeling, grid issues related to 
resolving near w'ake flow features such as the free shear 
layer persist 10 . 

The present experiments were performed in an 
effort to provide benchmark data for complex 
laminar/turbulent separated flows. Comparison of 
laminar heating prediction (via the LAURA code) to 
experimental measurement along the MSRO afterbody 
near impingement is showrn in Fig. 23, for M & =6 CF 4 
(Re 2D = 3,946). Measurement and prediction correspond 
to an open base. As discussed previously, the abscissa 
of this figure refers to the running length along the 
cylinder non-dimensional ized by the forebody diameter, 
originating at the aft comer and running forward 
towards the aerobrake base. The laminar experimental 
results indicate an impingement heating peak of 
h/h re j=0.4l at x/D=0.02 in close proximity to the 
afterbody comer. Numerical results obtained on both 
grids indicate a local heating spike at the afterbody 
comer (x/D=0) where there is a rapid thinning of the 
boundary layer (the spatial resolution of the 
thermography technique was insufficient to resolve a 
high gradient heating peak at this location), followed 
by a broad heating plateau (0.02<x/D<0. 1 ) interpreted 
as shear layer reattachment. As indicated in Ref. 39. 
the original grid used in this numerical simulation was 
adapted from a mesh utilized for early flight predictions. 
This modified grid had finer resolution in the near wake 
relative to that used for the flight cases but no attempt 
was made to align the grid w ith the shear layer. Some 
degree of numerical diffusion of the shear layer was 
expected. The numerical heating spike at x/D=0.24 is 
grid induced (see Ref. 39 for details). 

In an attempt to assess the sensitivity of wake 
closure characteristics and shear layer diffusion to grid 
alignment, a second calculation was made with the 
mesh realigned with wake streamlines. As anticipated, 
the realigned grid reduces the level of numerical shear 
layer diffusion and produces a sharper heating spike 
(h/h re p0.33) at impingement. The location of this 
heating maximum (x/D=0.12) occurs closer to the 
aerobrake base and further away from the experimental 
peak. 

Recent calculation of the Knudsen number near 
flow separation at the MSRO shoulder (not shown) 
suggest that, locally, application of slip boundary' 
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conditions at the wall near flow’ separation may be 
required. If practical, a 3-D coupled CFD-DSMC 
solution method could also prove beneficial in 
resolving differences between measurement and 
prediction. The present numerical predictions were 
made in the absence of a support sting. Computations 
with the support hardware modeled w'ould eliminate the 
uncertainty posed by this difference. In addition, 
comparison of the turbulent 1^=6 air data to numerical 
prediction using a two equation turbulence model are 
being considered. 

Concluding Remarks 

The Langley 20-Inch Mach 6 Air and CF 4 
Tunnels were used to assess the afterbody heating 
characteristics of a proposed Mars Sample Return 
Orbiter (MSRO). The Mach 6 CF 4 tunnel utilizes a 
heavy gas ( tetrafl uoromethane-CF 4 ) in order to generate 
a higher value of normal shock density ratio (than can 
be achieved in air tunnels) which is more appropriate to 
simulation of hypervelocity flight through planetary or 
earth atmospheres. The primary experimental 
objectives were to determine if shear layer impingement 
would occur on the orbiter afterbody at incidence angles 
appropriate to an aerocapture maneuver, and it present, 
determine the location and magnitude of the heating 
peak for comparison of CFD predictions to 
measurement. Global heat transfer mappings, surface 
streamline patterns, and shock shapes were obtained for 
post-normal shock Reynolds numbers (based on 
forebody diameter) ranging from 1,420 to 415,000, 
angles of attack ranging from -5 to 10 degrees at 0, 3. 
and 6 degree sideslip, and normal-shock density ratios 
of 5 and 1 2 . Laminar, transitional, and turbulent shear 
layer impingement on the cylindrical afterbody resulted 
in a localized heating maximum that ranged from 40 to 
75 percent of the reference forebody stagnation point 
heating. 

The present MSRO laminar data is consistent 
with a previously developed correlation for blunt body 
wake shear layer turning and has suggested the general 
applicability of the correlation towards initial design 
estimates defining the usable volume tor payload 
placement behind an aerobrake. If wake shear layer 
transition in flight is anticipated (ie. via forebody 
ablation or surface roughness) a larger shear layer 
turning angle should be expected. 
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Table 1. 

Nominal flow 

conditions 

ReJft(xlO h ) 

^ e 2.D 


Test gas 

P,.i(P*i) 

0.69 

36,700 

5.8 

air 

34 

1.23 

65,000 

5.9 

air 

64 

2.34 

123,000 

5.9 

air 

130 

3.32 

173,567 

5.9 

air 

186 

4.54 

235,967 

6.0 

air 

256 

6.19 

325.000 

6.0 

air 

372 

8.00 

415,000 

6.0 

air 

480 

0.01 

1.420 

6.2 

cf 4 

61 

0.02 

2,510 

6.0 

cf 4 

82 

0.03 

3,946 

6.0 

cf 4 

118 

0.13 

15,727 

5.9 

cf 4 

390 

0.23 

30,200 

5.8 

cf 4 

747 

0.36 

44,000 

5.9 

cf 4 

1073 

0.41 

52,442 

5.9 

cf 4 

1314 


T„(°R) 

862 

879 

902 

900 

900 

931 

926 

1266 

1255 

1238 

1210 

1281 

1248 

1277 


5.23 

5.24 

5.25 

5.26 

5.27 
5.27 
5.27 
12.1 
1 1.9 
1 1.9 

1 1.7 
1 1.8 
1 1.7 

11.8 
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Figure 7. Typical surface streamlines on MSRO 
afterbody in the presence of shear layer impingement. 
M*=6, air. Re, n =4 15,000, a= -4 degrees. 
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Figure 8. Typical surface heating on MSRO 
afterbody in the presence of shear layer impingement. 
M m = 6, air, Re :i> =4 15,000, a= -4 degrees. 



Figure 9. Effect of Reynolds number on normalized 
impingement zone heating along MSRO afterbody. 
M^=6, CF 4 , p 2 /p* = 12, a= -4 degrees. 



Figure 10. Effect of Reynolds number on normalized 
impingement zone heating along MSRO afterbody. 
M*=6, air. p : /p.„=5, ot= -4 degrees. 
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Figure 11. Effect of Reynolds number on normalized 
impingement heating peak on MSRO afterbody. 
M^=6, air and CF 4 , a= -4 degrees. 
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Figure 12. Comparison of present results to a unified 
wake transition correlation (Ref. 38). M,<=6, CF 4 , 

a= -4 degrees. 
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Figure 13. Effect of support sting diameter on 
normalized MSRO transitional impingement zone 
heating. M^=6, air, Re :n =36.700. a= -4 degrees. 



x/D 


Figure 14. Effect of open and closed MSRO base 
plane on normalized laminar impingement zone 
heating. M=6, CF 4 , Re :n =3,946, a= -4 degrees. 
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Figure 16. Effect of open and closed MSRO base 
plane on normalized turbulent impingement zone 
heating. M^=6, air, Re :{) =4 15,000, (X= -4 degrees. 



Figure 17. Surface streamlines on closed MSRO 
base plane in the presence of turbulent shear layer 
impingement. M„=6, air. Re 2n =4 1 5,000, 
a~ -4 degrees. 
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Figure 15. Effect of open and closed MSRO base 
plane on normalized transitional impingement zone 
heating. Mj=6. air, Re :D =52,442, a= -4 degrees. 


Figure 18. MSRO shear layer turning angle as a 
function of effective angle of attack. M c .=6, air. 
Re 2D =36.700 to 415,000. 
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